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We experimentally observe a counterpropagating dipole-mode vector soliton in a photorefractive SBN:60Ce 
crystal. We investigate the transient formation dynamics and show that the formation process differs 
significantly from the copropagating geometry. The experimental results are compared with fully anisotropic 
numerical simulations which show good qualitative agreement. @ 2008 Optical Society of America 
OCIS codes: 190.5530,190.5330 
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Spatial optical solitons in a counterpropagating (CP) 
geometry have been proposed some time ago.^ It has 
not been until recently however, that they have drawn 
wider attention^^ and were investigated theoretically in 
saturable Kerr and local photorefractive media in one 
transverse dimension (ID). In addition, Cohen et al. ex- 
perimentally demonstrated the existence of stable CP 
solitons with narrow stripe beams. ^ 

In a previous publication, we numerically predicted 
the existence of stable two-dimensional higher-order 
counterpropagating vector solitons in saturable Kerr me- 
dia,^ commonly used as an approximation for the pho- 
torefractive nonlinearity. In the present paper we will ex- 
perimentally and numerically demonstrate the existence 
of a counterpropagating dipole-mode vector soliton in 
a photorefractive SBN:60Ce crystal and investigate the 
transient formation dynamics. 

A dipole-mode vector soliton consists of two mutually 
incoherent beams: an optical dipole and a fundamental- 
mode(FM) beam. The individually propagating dipole 
does not form a spatial soliton due to repulsion of the 
dipole components.^' ^ However, if an FM beam which 
is incoherent to the dipole is launched in between the 
dipole components, they will be "trapped" via incoher- 
ent attraction.^' ^ 

The counterpropagating geometry poses some addi- 
tional challenges to the experimental setup. A soliton 
propagating inside a photorefractive crystal will be dis- 
placed in the positive direction of the crystal c-axis. This 
effect is well-known as the so-called "beam-bending" - 
effect caused by charge-carrier diffusion. This effect has 
been extensively investigated in the case of single soli- 
ton formation. For the case of interaction copropating 
solitons or copropagating vector solitons, the effect does 
not change the fundamental interaction behaviour as the 
beams are displaced by the same amount. Two coun- 
terpropagating solitons however are bend in the same di- 
rection. Thus when launched in a head-on configuration 



they will cross inside the crystal, resulting in a change of 
the soliton interaction length. Therefore, beam-bending 
has a significant effect on the transient formation dynam- 
ics as will become apparent later. Although at a first, 
rough glance the counterpropagating geometry might be 
very similar to configurations for e.g. double phasecon- 
jugate (DPC) mirrors, the actual experiments are quite 
distinct. Both interacting beams are mutually incoher- 
ent, thereby not interfering in the material. The lack 
of an interference grating is obvious due to the insensi- 
tivity of the beams to slight angular changes once the 
solitons are formed. This effect can directly be seen 
when observing the dynamics of the soliton formation. 
In contrast, the existence of a grating would result in 
a highly selective Bragg-condition. Such a behaviour is 
typical for double phase conjugation configurations, but 
does not appear in our case. Moreover, the scales are 
much smaller than for DPC experiments (input beams 
are in the range of 20 micrometers, interaction angles are 
smaller than 1°), and effects that are typically associated 
with four- wave mixing and phase-conjugation do not ap- 
pear on that scale. Nevertheless, in all our experiments 
we checked thoroughly that the effects we are observing 
are indeed caused by counterpropagating solitons. 

The experimental setup is depicted in Fig. Q A 
frequency-doubled Nd:YAG-laser beam is divided into 
two arms. The beams are rendered mutually incoher- 
ent by a means of a piezo-driven mirror (PM) oscillat- 
ing much faster than the response time of the crystal 
and are then focused onto the two faces of the pho- 
torefractive SBN-crystal via two beamsplitters. Due to 
the mutual incoherence of the beams, they do not inter- 
fere in the material nor produce an interference grating, 
thereby only interacting via their common phase modu- 
lation and induced refractive index channel. The crystal 
is biased with an external DC field along the x-axis par- 
allel to the crystallographic c-axis. To take advantage of 
the large vss electro-optic coefficient of SBN, the beams 
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Fig. 1. Experimental setup for the realization of a coun- 
terpropagating dipole-mode vector soliton. (PH: pinhole, 
F: lens, A/2: waveplate, (P)BS: (polarizing) beamsplit- 
ter, GP: glass-plate, M: mirror, PM: piezo-driven mirror) 



are polarized along the same direction. Background il- 
lumination is realized with a white light source, which 
allows adjustment of the dark-intensity. To be able to 
observe the two crystal faces simultaneously, both are 
imaged onto the same CCD-camera. The optical dipole 
is generated by entering a tilted glass-plate into one of 
the arms, a well-known technique for the creation of op- 
tical dipoles, and is oriented perpendicular to the c-axis 
to minimize effects of the photorefractive anisotropy.^^ 
In order to compensate for the above-mentioned beam- 
bending, the beams are launched into the crystal in such 
a way that the input face of one beam coincides with the 
exit face of the counterpropagating beam propagating as 
an individual soliton. This results in a small angle (< 1°) 
between the two incident beams. The power of the FM 
beam is PFM-beam ^ 1-00 jiiW and the power of the 
dipole Pdipoie ^ 1-35 /j^W. The approximate diameter of 
the beams at the input face is (i 20 /im FWHM. The 
geometry of the SEN crystal is 23(a) x 5(6) x 5(c) mm 
and propagation is along the a-axis of the crystal for 
maximum propagation length, corresponding to approxi- 
mately 4 diffraction lengths. The external field was set to 
Eext « 1.2 kV/cm and the background-illumination was 
adjusted in order to allow stable self-focusing and soliton 
formation of the individually propagating FM beam. 

The formation of the dipole-mode vector soliton is 
shown in Fig. [21 The images are taken at the exit face 
of the dipole beam, the less intense spot in the mid- 
dle of the images is the reflection of the FM input at 
the crystal face. It serves as a reference and does not 
influence the soliton formation process. When the non- 
linearity is turned on, the dipole initially self-focuses to 
the left of the FM beam (Image (b)). This is due to 
the initial angle between the beams. The dipole is then 
strongly deflected to the right, as beam-bending and at- 
traction from the FM soliton act in the same direction. 
The dipole splits up and two spots of low intensity with 
a much smaller separation appear on the left side of the 
FM beam (c). As time elapses, this trapped part of the 



Fig. 2. Images of the dipole exit face of the crystal(beam 
from FM input face to camera is blocked). The spot in 
the middle is the reflection of the FM input serving as 
a reference, (a) linear propagation. Soliton formation at 
times t = 25 s (b), t = 128 s (c) and t = 500 s{d) af- 
ter the nonlinearity has been turned on. (e) stable state 
after t ^ 35 min. (f) individually propagating dipole ap- 
proximately t = 5 min after the FM beam was blocked. 



dipole draws intensity from the non-trapped part and 
the non-trapped part is attracted horizontally toward 
the FM input (d). After approximately 35 minutes, the 
system reaches equilibrium, with almost all power of the 
dipole trapped (e). When the FM beam is blocked, the 
dipole components quickly repel and move to the right 
due to beam-bending. Image (f) shows the dipole after 
propagating individually for approximately 5 minutes. 
The separation of the dipole components of the trapped 
and non-trapped part in image (e) was dtrapped ^ 38 fim 
and dnon-trapped ^ 87 /im respectively, compared to 
dindividuai ~ 57 /im for the individually propagting 
dipole. 

The transient formation dynamics differ significantly 
from the copropagating case, which makes them partic- 
ularly interesting. The most striking fact is the splitting 
of the dipole into two parts, which is not observed in 
the copropagating geometry. Although at a quick glance 
such an effect might be attributed to phenomena as con- 
ical scattering, a thourough analysis shows immediately 
that those effects can be excluded: the two interacting 
solitons are mutually incoherent, and thereby not able to 
interfere with each other. Moreover, the angles at which 
these new beams appear are in the range of one degree, 
much smaller than the scattering cone which should be 
well above 10°. Instead, the splitting can be explained 
by the angle between the beams which is introduced to 
compensate effects of beam-bending. Soliton formation 
in photorefractive crystals is a two-step process, in which 
the beam self-focuses initially and is then displaced by 
beam-bending.^^ After self-focusing of the counterprop- 
agating beams self-focus but before beam-bending takes 
place, the beams only cross at one point inside the crystal 
due to the angular adjustment. At this point the colliding 
beams will interact strongly, which causes the observed 
splitting. The non-trapped part of the dipole is then de- 
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Fig. 3. Images of the dipole input face at steady state 
^ 35 min after the nonhnearity has been turned on. (a) 
The output of the FM beam is displaced to the left of 
the dipole input, (b) Output of the FM when the dipole 
is shortly blocked, a small fraction of the FM beam has 
coupled into the dipole waveguide. 
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Fig. 4. Numerical simulations of the dipole-mode vector 
soliton. The images show the exit face of the dipole at 
(a) t = 100 r and {h)t = 280 r after the start of the 
simulation. The location of the input beams is indicated 
by the dashed lines. 



fleeted by the combined forces of beam-bending and in- 
coherent attraction of the dipole, causing the dipole to 
"overshoot". This explains why the non-trapped dipole 
part is initially displaced further then the individually 
propagating dipole, which can be seen by comparing im- 
ages (c) and (d) with image (f) in Fig.[2j The second re- 
markable feature of the counterpropagating dipole is the 
long timescale of the formation process. It differs signifi- 
cantly from copropagating solitons, with formation times 
in the range of 20 — 50 seconds or even shorter for higher 
intensities,^^ compared to around 30 minutes for the 
counterpropagating dipole-mode vector soliton to reach 
a stable state for similar parameters. As can be seen from 
Fig. [2fc)-(e), the trapped dipole is not aligned with the 
input of the counterpropagating beam, but is displaced 
slightly to the left. The output of the FM beam, which 
can be seen in Fig. IHl shows a similar picture. The FM 
soliton is displaced to the left of the dipole input. This 
is in agreement with earlier simulations^ which predict a 
threshold propagation length where the symmetry along 
the dipole axis is broken. Thus our experimental results 
are above this threshold length. Additionally we observe 
that the the FM beam splits as well. Fig. ISfb) shows 
the exit face of the FM beam when the dipole is shortly 
blocked. A fraction of the beam has split from the main 
part and has coupled into the dipole waveguide. 

To verify our experimental results, we carried out sim- 
ulations in the anisotropic, nonlocal model. The two 
counterpropagating beams are denoted by their slowly 



varying envelopes F and B. The total optical field can 
then be written as F ex.p{ikz) -\- B ex.p{—ikz)^ with k be- 
ing the wavevector in the undisturbed crystal. Since we 
are considering F and B to be mutually incoherent, the 
total intensity is just given by / = |Fp + The 
Kukhtarev model of the photorefractive crystal leads un- 
der simplifying but well justified assumptions to the fol- 
lowing equation 

-^dtAcj) + A(/) + V ln(l + /) V0 = 

^exAln(l+/) 
-^f^[A ln(l + /) + (V ln(l + (1) 

for the potential (f) of the electrical screening field, which 
causes the nonlinear change of the refractive index via 
the Pockels effect. E^xt denotes the externally applied 
electrical field, = —kBT/e is the diffusive coupling 
strength, r is the effective time constant of the crystal, 
and / is measured in units of the background intensity. 

The propagation of the beams is in paraxial approxi- 
mation described by the set of equations 

dzF - '-VlF + '-i{E,,tF - dx^F) = (2a) 
-d,B - '-VlB + '--i{E,,tB - dx(^B) = (2b) 

with the nonlinear coupling constant 7 = /c^Xongreff, 
which contains the refractive index of the unperturbed 
crystal no, the effective element of the electro-optic ten- 
sor Teff, and the transverse scaling constant xq. de- 
notes the transverse Laplacian. The propagation (2:-) axis 
is scaled to the diffraction length Ljj = kx^. The values 
used in the simulations are xq = 10 /im (typical beam 
widths are about 1 . . .2xo), rgff = 280pm/V, no = 2.35, 
Eext = 2.5kV/cm, and k = 21:710/ \ with A = 532 nm. 
The total propagation distance was chosen to be AL^ 
(11.1 mm). 

The results of the simulations are shown in Fig.^ Im- 
age (a) and (b) show the exit face of the dipole beam 
at times t = 100 r and t = 280 r after the start of 
the simulation. Although the numerical and experimen- 
tal results differ in detail, the main features of the in- 
teraction are very similar. The splitting of the dipole 
into a trapped and non-trapped part can clearly be ob- 
served in Fig.^a), although it is not as pronounced as in 
the experiment. As time elapses the trapped part grows 
in intensity drawing energy from the non-trapped part 
(Fig. 2|Jb)). The second main feature, the exceptionally 
long timescale is reproduced as well, with the time con- 
stant r corresponding to approximately r ^ 15 s. 

To summarize, we have shown the existence of a sta- 
ble counterpropagating dipole-mode vector soliton in a 
photorefractive SBN crystal. The vector soliton differs 
considerably in many aspects from its counterpart in co- 
propagating geometry. The timescale of the transient dy- 
namics are significantly larger. During the formation the 
beams split up, and a trapped and non-trapped part of 
the dipole can be observed. 
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